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Abstract. We compare the observed spatial offsets of short gamma-ray bursts from their 
host galaxies with their predicted distributions, assuming that they originate in double neu- 
tron star binaries that form from field stars. We find that, for the majority of bursts, this 
model is sufficient to explain the observed offsets, although there is a trend towards larger 
offsets than predicted. One burst, GRB 060502B, has an offset that is clearly anomalous. We 
discuss possible reasons for the large offsets, including host galaxy misidentification, and 
suggest that some of the largest-offset bursts may originate in the merger of double neutron 
star binaries that form dynamically in the cores of globular clusters. 

Key words, binaries: general - gamma-ray burst: general - neutron stars: general - black 
holes: general 



1. Introduction 

The advent of the Swift satellite has, for the 
first time, allowed us to locate the positions 
on the sky of short-duration gamma-ray bursts 
(SGRBs) and hence determine their host galax- 
ies. The sample of located short bursts that has 
built up over the last six years of Swift oper- 
ations covers a range of galaxy types, includ- 
ing elliptical galaxies in which no stars are cur- 
rently being formed. The bursts are observed to 
occur with a wide range of spatial offsets from 
their inferred host galaxies, in some cases oc- 
curring well outside the host. 

These observations broadly support a pic- 
ture in which the bursts are powered by the 
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formation of a stellar-mass black hole during 
the merger of a binary containing either two 
neutron stars (NS-NS) or a black hole and 
a neutron star (BH-NS). These mergers are 
driven by inspiral caused by the emission of 
gravitational radiation; this leads naturally to a 
very wide range of merger timescales. Hence 
SGRBs can occur in non-star-forming hosts. 
Meanwhile, the kicks expected to be present at 
the formation of the neutron stars will impart 
a natal velocity to the binaries, which offers 
a natural qualitative explanation of the offsets 
(e.g. Bloom et al. 1999; Fryer et al. 1999). 

Utilising the hosts identified by Swift, we 
build on the work of previous authors by con- 
sidering the bursts on a host-by-host basis. 
We construct offset distributions based on the 



Church et al.: SGRBs from globular clusters 



1 



GRB 


^off 




Vh 


fh 


Rz 




kpc 




km/s 


kpc 


kpc 


050509B 


63.7 


± 12.1 


664 


46.3 


21.0 


050709 


3.55 


±0.27 


110 


7.9 


1.8 


050724 


2.54 


±0.08 


532 


23.9 


4.0 


051221A 


1.53 


±0.31 


157 


15.7 


2.2 


060502B 


73 ± 


13 


505 


20.5 


10.5 


060801 


19.7 


± 14.0 


170 


18.2 


3.0 


061006 


1.44 


±0.29 


124 


9.9 


3.7 


061201 


33.9 


± 0.4 


121 


9.6 


1.8 


061210 


10.7 


± 6.9 


162 


16.5 


2.6 


061217 


55 ± 


20 


141 


12.8 


1.8 


070429B 


4.7 ± 


4.7 


149 


14.2 


2.1 


0707 14B 


3.08 


± 0.47 


111 


8.9 


0.94 


070724A 


4.76 


±0.06 


435 


13.1 


3.2 


070809 


19.61 


± 1.9 


110 


8.0 


0.92 


071227 


16.1 


± 0.2 


173 


18.8 


3.1 


080905A 


18.11 


± 0.42 


170 


18.1 


3.0 



Table 1. Properties of SGRBs in our sample 
and their host Galaxies. For more details of the 
population and a complete reference list see 
(Church etal. 2011). 



properties of the observed hosts. We use a sam- 
ple of 16 bursts, all of which have identified 
hosts with measured properties including mag- 
nitude and redshift. Our burst sample is pre- 
sented in Table 1 . 

2. Data and calculations 

In order to predict the offsets of short gamma- 
ray bursts we synthesised a large population of 
compact binaries using the rapid binary evo- 
lution code BSE (Hurley et al. 2002). We re- 
tained those binaries that evolved into NS- 
NS or BH-NS binaries, and computed their fi- 
nal 3D velocities, taking into account the ef- 
fects of supernova mass loss and natal kicks. 
We found that, in order to satisfactorily re- 
produce the orbital properties of Galactic NS- 
NS binaries, we required strong natal kicks, of 
the order of 100 km s . Hence we adopt the 
Arzoumanian et al. (2002) kick distribution. 
We present distributions of merger times and 
rest-frame velocities for our final sample of bi- 
naries in Figures 1 and 2. 

We model the observed SGRB sample on 
a burst-by-burst basis. For each host galaxy 
we produce a separate potential model, utilis- 
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Fig. 1. The merger time distribution of our synthe- 
sised binary population. NS-NS binaries are shown 
in as solid lines; BH-NS binaries as dashed lines. 
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Fig. 2. The velocity distribution of our synthesised 
binary population after the formation of a dou- 
ble compact binary. NS-NS binaries are shown in 
as solid lines; BH-NS binaries as dashed lines. 
Velocities are measured with respect to the initial 
rest frames of the binaries. 



ing the logarithmic potential of Thomas et al. 
(2009). The core radius and halo circular 
velocity Vh for each burst are obtained from 
the the fits of Kormendy & Freeman (2004) 
to SDSS data. Fits from the same source were 
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Fig. 3. The observed offset from the host centre for 
each of the bursts in our sample. To allow the con- 
sistent comparison of different hosts we have plotted 
the difference between the observed offset and the 
median predicted offset, in units of the standard de- 
viation of its host's predicted offset distribution. For 
example, a burst which lay at the 98th percentile in 
the cumulative offset distribution of its host would 
be plotted at x = 2. If the predicted distributions 
match the observed values then the plot should be 
symmetric around x = 0. For each burst the range 
plotted is one standard deviation in observed offset. 
The lower bar in each case is for NS-NS progeni- 
tors, the upper for NS-BH progenitors. 



also used for half-light radii where measured 
values were not available in the literature. The 
properties of our haloes are given in Table 1. 
We form stars in an exponential disc with ra- 
dial scale length equal to the half-light radius 
of the host. We place our compact binaries 
in these discs, endowing them with the calcu- 
lated velocities in isotropically distributed di- 
rections. We integrate the motion of each bi- 
nary's centre of mass in the potential field of 
the host for the merger time of the binary, 
recording the position in the host galaxy where 
the merger occurs. Finally we project the loca- 
tion onto the host using a random viewing an- 
gle. This process is repeated over a large num- 
ber of realisations to build up an offset distri- 
bution for each host. 

The position of each observed burst within 
the offset distribution predicted for its host 
is plotted in Figure 3. The majority of the 
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Fig. 4. Lower panel: the cumulative offset proba- 
bility for bursts around the host of GRB 060502B 
assuming that the burst comes from a merging field 
binary. Solid lines show the distributions for bursts 
from NS-NS binaries; dashed lines show bursts 
from BH-NS binaries. The error bar, placed at ar- 
bitrary height, shows the 1-cr error on the mea- 
sured offset. The dotted line shows the 3-<x exclu- 
sion level. Upper panel: ditto, except for NS-NS bi- 
naries that form dynamically in a globular cluster 
system similar to that around M87. 



bursts are reproduced relatively well by our 
treatment, although it is evident that the syn- 
thesised distributions are systematically under- 
predicting the host offsets; i.e. the distribu- 
tion is not centred around x = 0. In particu- 
lar one burst, GRB 060502B, is at an anoma- 
lously large offset. The cumulative distribution 
of predicted offsets for this burst is shown in 
Figure 4. The possibility of a burst occurring 
at such a large offset around this host galaxy is 
excluded at the three-sigma level, although the 
errors on the burst position are rather large as 
no optical afterglow was detected. 

3. Discussion 

There are several possible explanations for the 
failure of the field binary scenario to reproduce 
the large offset of GRB 060502B. Firstly, it 
is possible that the host galaxy has been mis- 
identified. We follow Bloom et al. (2007) in 
selecting the giant elliptical considered here 
as the host with a large offset. The probabil- 
ity of a chance co-incidence with this galaxy 
is only a few per cent; on the other hand there 
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are smaller, fainter galaxies within the XRT er- 
ror circle (Berger et al. 2007). This is an innate 
problem in measuring offsets with hosts se- 
lected through spatial co-incidence on the sky. 

A second possibility is that, for bursts 
around massive elliptical galaxies, the evolu- 
tion of the hosts has driven the binaries out 
to larger offsets. Zemp et al. (2009) compute 
the distribution of coalescing compact binaries 
in cosmologically evolving dark matter halos 
and show that the result is to increase the ob- 
served offsets, with the inferred host galaxy not 
always being the galaxy in which the binary 
formed. This is a entirely possible origin for 
the progenitor of GRB 060502B. 

The third possibility we consider is that 
some fraction of bursts originate in com- 
pact binaries that have formed dynamically 
within the globular cluster systems of their 
host galaxies. Within the cores of globular 
clusters the number densities are high enough 
that stars undergo dynamical encounters with 
one another. Neutron stars, which segregate to 
the centre of the clusters, can exchange into 
pre-existing binaries to form NS-NS binaries 
(Davies 1995). Those binaries then could yield 
SGRBs. Grindlay et al. (2006) extrapolate 
from M15-C, the single NS-NS binary that is 
unambiguously dynamically formed, to derive 
a merger rate of R gc = 0.3ir(10 9 M o )- 1 f H ^ bble , 
where T is the number of globular clusters 
per 10 9 M© of galactic stellar mass. They cau- 
tion, however, that simple scattering calcula- 
tions suggest that this rate is an underestimate 
by a factor of at least ten. Our population syn- 
thesis implies a field NS-NS merger rate of 
Afield = \5M\0 9 M e )- 1 

bubble' wri ere fi, is the 
field binary fraction. This suggests that the ra- 
tio of rates of field and globular cluster merg- 
ers is 0.0277/t,. The nearest ellipti- 
cal galaxy of similar size to the putative host 
of GRB 060502B is M87, which has roughly 
14 000 globular clusters (Harris 2009) in a very 
extended distribution. Its specific density of 
globular clusters is T 8 (Brodie & Strader 
2006). Taking into account the underestimate 
mentioned above and the uncertainty in binary 
population synthesis calculations this is consis- 
tent with a similar rate of burst production in 
the field and in globular clusters for massive el- 



liptical galaxies. In the upper panel of Figure 4 
we plot the cumulative offset distribution given 
by assuming that GRB 060502B occurred in a 
globular cluster in a system identical to that of 
M87. This gives a much better fit to the offset 
than a field population origin. 
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